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Abstract 

Bounds on the exponential decay of generalized eigenfunctions of bounded and un- 
bounded selfadjoint Jacobi matrices in ^^(N) are established. Two cases are considered 
separately and lead to different results: (i) the case in which the spectral parameter lies 
in a general gap of the spectrum of the Jacobi matrix and (ii) the case of a lower semi- 
bounded Jacobi matrix with values of the spectral parameter below the spectrum. It is 
demonstrated by examples that both results are sharp. 

We apply these results to obtain a "many barriers-type" criterion for the existence 
of square-summable generalized eigenfunctions of an unbounded Jacobi matrix at almost 
every value of the spectral parameter in suitable open sets. In particular, this leads to 
examples of unbounded Jacobi matrices with a spectral mobility edge, i.e. a transition 
from purely absolutely continuous spectrum to dense pure point spectrum. 

1 Introduction 

One of the central tools in the spectral theory of differential and finite difference operators, in 
particular Schrodinger operators and their discrete counterparts, are results on the asymptotic 
behavior of generalized eigenfunctions. Here we are concerned with establishing such results 
for unbounded Jacobi matrices and relating them to spectral properties of the associated self- 
adjoint operators. 

In the first part of the paper we consider general unbounded self-adjoint Jacobi matrices 
J on £^(N) and assume that A G M lies in a spectral gap of J. We will use a "discrete" and 
rather simple version of a technique introduced in [1] to prove upper bounds on the exponential 
decay of generalized eigenfunctions of J to A. The decay bound for eigenfunctions of Schrodinger 
operators found in improved on longstanding bounds obtained through the so-called Combes- 
Thomas method |1] . Our results on exponential decay are stated in Section [2] and proven in 
Section [31 

While these results are quite general and their proof, due to the discrete one-dimensional 
setting, quite elementary, the obtained bounds are remarkably sharp in several respects. This 
will be understood in Section HI where we will consider two concrete classes of unbounded 
Jacobi matrices for which the exact asymptotics of generalized eigenfunctions can be obtained. 

Combes-Thomas type estimates like the ones proven here are a frequently used tool in the 
spectral analysis of differential and finite difference operators. Some new applications are given 
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in Sections [5] and [6] below. We also mention a recent paper by J. Breuer [2], where a similar but 
weaker estimate for the matrix elements of (J — z)^^ was found and applied to study spectral 
and dynamical properties of random Jacobi matrices. 

Our main application of the results on exponential decay is a criterion for the existence of 
£^-solutions of Ju = Xu at (Lebesgue) almost every energy A in an interval I. This describes 
types of Jacobi matrices J which coincide on infinitely many growing intervals with a Jacobi 
matrix Jq which has J as a spectral gap. Physically, these intervals can be thought of as a 
series of barriers preventing quantum mechanical transport under the time evolution for the 
hamiltonian J. Consequences are the absence of absolutely continuous spectrum for J in / 
and, after adding generic rank-one perturbations, pure point spectrum in /. For discrete and 
continuous one-dimensional Schrodinger operators results of this kind were first obtained by 
Kirsch, Molchanov and Pastur in [151 HB]- Subsequently, their ideas have been extended to 
multi-dimensional Schrodinger operators [TJ EH], as well as to finite difference operators on 
strips [T7] . 

We state and prove a general barriers-type criterion for unbounded Jacobi matrices in 
Section O In the final Section [HI we study a specific class of unbounded Jacobi matrices, which 
give rise to a mobility edge. Their spectrum covers the entire real line, it is purely absolutely 
continuous outside an interval [— c, c] for some c > 0, and dense pure point in (— c, c). The latter 
will follow from our criterion in Section [5] together with a Weyl-sequence argument, while the 
claim on absolute continuity will be a consequence of the general results on the asymptotics of 
solutions of difference equations in [S]. 

2 Exponential decay bounds for the resolvent 

Let (•, •) be the inner product in £^(N) and denote by e„ the canonical basis in £^(N). Before 
we proceed further recall a general result from [TB|, going back to [S] and [T^, which says that 
for any bounded operator A in £^(N) with the band matrix aij = {Aej, Cj) having the bounded 
inverse A~^, the entries bij of A"^ satisfy the estimates: 



for some C > and r G (0, 1). In the proof of this estimate the boundedness of A is used in an 
essential way. 

Consider a Jacobi matrix J defined by the difference relations 



hil < Cr 





and boundary condition u{0) = 0, or, in equivalent matrix notation. 



/ ^1 
Ai 



q2 A2 
A2 q-s 



\ 



J 




\ 
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Here, the "discrete potential" g„ is real and the "weights" A„ are positive for n G N. For later 
use we also adopt the convention Aq = 1. Sometimes it is also convenient to represent J in 
the form J = SA + AS* + Q, where A (respectively Q) are the diagonal operators given by A^ 
(respectively g^) in the canonical basis e„, n G N , and S is the unilateral shift S'e„ = e„+i. 
Under the assumption 

oo ^ 

5^- = 00 (2.3) 

n=l 

this defines a unique self-adjoint operator on ^^(N), e.g. [3], which will also be denoted by J. 

The main result of the first part of this paper is an upper bound for the matrix elements of 
the resolvent of J: 



Theorem 2.1 Let J* = J and assume that liminf„,_»oo A„ > 0. 

(a) Let (r, s) be a finite gap in the spectrum of J. Then there exists a constant f] > such 
that for all A G (r, s) and all n G N, 



n— 1 



|((J - A)-iei, en)\ < 4max {(A - r)-\ {s - A)-^} exp ( -^^{X - r){s - X)J2 y ] ■ (2-4) 



k=l 



(b) Suppose that J is bounded from below and denote d = inf cr{J). There exists a constant 
Tj > such that for all X G (—00, d) and all n eN, 



n—l 



|((J-A)-^ei,e„)| <4(t/-A)-iexp \^-^y^d^Yl—j ■ (2.5) 
The corresponding result holds if J is bounded above and X G (sup (t( J), 00). 



In the case lim„^oo A„ = +00 the above estimates can be slightly improved. This is the 
content of the next two results. 



Theorem 2.2 Let J* = J. Suppose that lim^^oo A^ = +00 and let (r, s) be a gap in the 

spectrum of J. Then for arbitrary e G (0, 1/2) there exists N = N{e) such that 

- ^)-'^-»«>l S (-(i - .)V(A-0(-A)g±) , (2.6) 

for all X G (r, s) and for all n > N. 

For a Jacobi operator J which is bounded from below and A below the bottom of the 
spectrum the estimates given in the above results can be further improved. Indeed, we have 
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Theorem 2.3 Suppose that J = J* associated to the weights Xn and the potential qn is bounded 
from below by d and lim„^oo A„ = +00. Fix e G (0, 1) and complex A such that 3?A < d. Then 
there exists N = N{e, A) such that 

|((J-A)-iei,e„)| < [{d-^X)e\-'exp ( -(1 - e)Vd^^J2 4r] ' (2.7) 
for n > N. 

In Section H] we will demonstrate by examples that the constants in the exponents on the 
right hand sides of (12.61) and (12. 7p are optimal in the sense that 1/2 — e and 1 — e, respectively, 
can not be replaced by 1/2 + e and 1 + e for any e > 0. 

The sequence v defined by v{n) = ((J — A)~^ei,e„) is a generalized eigenfunction of J, 
meaning that it solves (12. ip for n > 2 without satisfying the boundary condition at 0. As 
J is in the limit point case at +00, it is (up to constant multiples) the unique generalized 
eigenfunction which is square-summable. Thus we may understand (12. ip as a bound on the 
decay of generalized eigenfunctions. As A ^ cr( J) (the spectrum of J), f is not an eigenfunction, 
but the above results also provide bounds on eigenfunctions for perturbations of J: If J = J+A, 
where A is a finite Jacobi matrix, and if u = {u{n)) is an eigenfunction of J for an eigenvalue 
A G (r, s), then u satisfies the bound given by the right hand side of (12. 4p . This is obvious as 
u{n) for large n coincides with the unique £2-generalized eigenfunction of J. 

We point out two specific features of the exponent on the right of (12.40 : It describes the 
large n asymptotics as well as the asymptotics as A approaches the spectrum of J. The large n 
asymptotics, characterized through the sum ^1/Afe, generalizes the "Schrodinger case" A^ = 
1, where generalized eigenfunctions for A in a spectral gap decay exponentially in n. As A 
approaches a spectral edge of J, that is, either r or s, the constant in front of ^1/Afc is 
proportional to the square root of the distance of A to the spectrum. This improvement over 
the original Combes-Thomas method, which merely provides a term which is linear in the 
distance, is due to the new ideas introduced in [Ij. We will comment on this at the end of the 
proof of Theorem 12.11 in Section [31 

The proofs of Theorems 12.11 to 12.31 allow for several generalizations. As an example, we 
state the following generalization of Theorem 12.1( a). which will be used in our applications in 
Section O Similar generalizations could be formulated for our other results. 

Theorem 2.4 Under the conditions of Theorem \2.1\ there exists rj > such that for all A G 
(r, s), all real 6 with \6\ < |a/(A — r)(s — A) and all subsets A,B C'N with maxB < minA it 
holds that 

(miriA— 1 ^ 
-W{X-r){s-X) V - 

(2.8) 

Here xa and xs denote the multiplication operators with the characteristic functions of A 
and B and || ■ || the operator norm. Theorem 12.11 is a special case of Theorem 12. 4[ where the 
energy is real {6 = 0), A = {n} and B = {1}. 



Unbounded Jacobi matrices 



5 



3 Proofs of Theorems EH, [2l2] and [2I3 



3.1 Proofs of Theorem D and [2741 

Let p{n) := X]fc=i V-^fc ^^^d, for 7 > to be specified later, 
operator in Then a calculation shows that 



J = A(7) 



/ ai 

61 a2 

bo 



-7P 



multiplication 



(3.1) 



is a non-symmetric Jacobi matrix with entries 



1 



bn — A„ 



n eN . 



(3.2) 



0(n) 7' ■•"^0(^ + 1; 

To determine bounds on the real and imaginary part of A we verify by Taylor expansion that 



a„ + bn 



(3.3) 



and 



We conclude that 



( 

ax + 61 



-27 + 0(7VAn') . 



(3.4) 



ai + 61 

02 + &2 

02 + &2 



r 
2 



/ Ar^ 

' A,-^ 



Ar 



A,^ 



and 



ImA(7) = i 



/ a\ — b\ 

b\ — a\ — 62 

62-02 

V 



/ 1 

-1 1 



«7 



-1 



V 



(3.5) 



-0(f ) , (3.6) 



where An > C > uniformly in n was used and error terms refer to norm bounds. In fact, all 
we will use below are the norm bounds ||ReA(7)|| < Ci7^ and ||Imyl(7)|| < C27 with Ci,C2 
depending only on J. 

The following lemma provides the operator theoretic fact behind the improvement of the 
Combes- Thomas method found in [1]. Extracting from the argument in [1], this lemma is stated 
with proof in [19], p. 60]. 
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Lemma 3.1 Let T be invertible self-adjoint with = dist(0, cr{T) fl (0, oo)) and 

c/_ = dist(0,or(T)n(-oo,0)). Let S be self- adjoint, \\S\\ < 1. Then for /3 e R, \(3\ < l^d+-d-, 

the operator T + if3S is invertible, with 

||(T + 2/35)-i <2-max|i-,i-| . (3.7) 

Remark 3.2 If a(T) C M+ (resp. o"(T) C M_ ) then we put d_ = oo ( resp. dj^ = oo ) . 

We apply Lemma [311] to the operator J — A + ^(7) = T + i[3S, with self-adjoint operators 
T = J- A + Rev4(7) and S = ImA(7)/||Im A(7)||, and p = ||ImA(7)||. As ||ReA(7)|| < C172, 
the operator T has a spectral gap (r — A + Ci7^, s — X — Ci7^) = (— (i_|_). 

For a finite interval (r, s) we choose 7 = //a/ (A — r)(s — A), where 

"-"■'■'{i^^- V2C.|.-r) }- 

Then d_ > i(A - r), d+ > i(s - A) and \(3\ < i^(A-r)(s-A) < ^^/d+d^. Applying 
Lemma [3. II we know that J — A + ^(7) is invertible and thus, using (13. 7p . 

||0-i(J-A)-V|| = ||(J-A + A(7))-^|| 

1 1 



< 2 max 



d. ' (i_ 



We note that |((J — A)~"^ei,e„)| = ||e„(J — A)""'^ei||, where on the right we think of e„ and ci 
as multiplication operators (and of || ■ || as the operator norm). The proof of Theorem 12.1( a) is 
thereby completed through the estimate 

||e„(J-A)-V|| < ||e„0||||0-i(J- A)-V||||0"'ei|| 

< e-^^(") ■4max|^-,-^l ■ 1 . (3.10) 
l^s — A A — rj 

Part (b) follows by simple modifications of the above argument. We are in the case d- = 00, 
meaning that there is no restriction on the size of 1/3 1 . One chooses 7 = ^{d-X)/2Ci. 

The proof of Theorem 12.41 also follows the above lines with only minor changes. The imag- 
inary part of the spectral parameter is included in 5* by choosing S = {61 + Imyl(7))/(|5| + 
||Imy4(7)||), and /5 = |5| + ||Im74(7)||. The constant rj is modified to be the minimum of 
and , The calculation (13.101) applies in the same way to ||xa(-^ — (A + i6))~^XB\\ to 

y/2Ci{s-r) 

give the bound (12.81) . 

Note that in the above argument crucial use was made of the fact that the real part of 
^(7) is smaller than its imaginary part, 0(7^) as compared to 0{'~f). This is exploited through 
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Lemma 13.11 which can be interpreted as saying that purely imaginary perturbations have a 
weaker effect on the invertibihty of self-adjoint operators than self-adjoint perturbations. In the 
"classical" Combes-Thomas method only the rougher bound ||y4(7)|| = 0(7) is used, yielding 
merely linear dependence of the exponent in fl2.4p on the distance of A to the spectrum. 

3.2 Proof of Theorem [2721 

We start with an improved version of Lemma 13. 1[ 

Lemma 3.3 Let T be invertible self-adjoint with positive = dist(0,cr(T) fl (0, 00)) and 
d_ = dist(0,cr(T)n(-oo,0)). Let S be self-adjoznt, \\S\\ < 1. Then for f3eR, \/3\ < .J d^ ■ 
the operator T + i(3S is invertible, with 

UT + ^(3S)-'\\ < [A+-(A^+/?2-rf+.rf_)V2]-i, A+ = l(rf_ + rf+). (3.11) 



Remark 3.4 If a{T) C M+ ( resp. cr{T) C M- ) then we take the limit of the right hand side 
of 1^3. as d^ tends to +00 (resp. as d+ tends to +00 ), and put no restriction on (3 . 

Proof: Let us start with some geometrical ideas which are behind the proof. Denote by K be 
the upper half-disc with the center at {d^ — dJ)/2 ( on the x-axis ) of radius A_|.. Then the 
length of the segment from (0, 0) to the point of the intersection of the upper circle with the 
y-axis equals ^Jd+ ■ d- . This explains why < \/d+ ■ d_ . 

Let P+ (respectively P_) be the spectral projection of T corresponding to the positive 
(respectively negative) part of o"(T). We have 

T = |r|(p+ - p_) = Ti(rf+p+ - rf_p_), 

where Ti := jT'K^ — ^)- Note that Ti commutes with P+ and P_ and Ti > /. 



Hence 



II (T + ipsy^W = \\T-^^\d+P+ - d^P^ + ipT~^'^ST~^'\^T^^'^\\ 
< \\{d+P+-d^P_+i(3Si)-% 

with SI = Si := T^^^^ST^^^\ \\Si\\ < 1. 

Therefore the proof is reduced to the case T = d^P^ — d_P_. Denote by A_ = {d^ — dJ)/2 
and Ji = P+ — P_. Then 

T + ij3S = A_I + A+ Ji + t/3S = A+[Ji + A;^(A_/ + i/3S)] (3.12) 

and 

II A_/ + if3Sf = \\All + f3^S*S\\ < +(3^ = Al + {/3^ - d_d+). (3.13) 
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Using (13.121) and fl3.13p we obtain 



\\{T + ipS)-'\\ < A;i[Ji + (A_/ + z/5^)A- 



n-ii 



This completes the proof of Lemma 13. 3[ 

Turning to the proof of Theorem 12.21 we fix e G (0, 1/2). Let 

n-l 

p(n) := Oforn < A^, p{n) := ^ l/A^ for n > N, 

k=N 



(3.14) 



with to be chosen below. For 7 > to be specified later, let = e~'^'' be a multiplication 
operator in £^(N). Then by repeating the calculations given in the proof of Theorem 12.11 we 
find 

/ ai \ 

bi a2 



/-I 



J0- J = A(7) 







(3.15) 



is a non-symmetric Jacobi matrix with entries 

(n + 1 



(p{n) 



- 1 



(p{n) 



(n + 1) 



We want to determine bounds on the real and imaginary part of ^4(7). Note that 

ReA(7) = SW + WS*, = Dm^(0,...,V'i(iV),^i(A^ + 1),...), 
ImA(7) = i{SZ-ZS*), Z = Diag{0,...,2P2{N),i:2{N+l),...), 



(3.16) 



(3.17) 
(3.18) 



where 

00 00 

Up) = 5^f1(2A:)!Af-i]-\ Mp) = -Y.i''-'[{2k - ly.xf -^' ■ 

k=l k=l 

We apply Lemma to the operator J — A + ^(7) = T + ipS, with self-adjoint operators 
T = J-A + Re A(7) and S = Im A(7)/||Im A(7)||, and P = ||ImA(7)||. Using (IXT7D and (IXT5D 

we have 



e(iV) = ||ReA(7)||, 5(iV) = ||ImA(7)||, 

where e(A^) := 2 sup^ t/'i (p) and S{N) := 2supp {—■ip2{p))- 

Let 7 = (1/2 — e)\/(A — r){s — A). Easy computations show that 



(3.19) 



6(iV)<Ci(J)(A-r)(.-A)(inf A,) 

p>N 



-1 



(3.20) 
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and 

S{N) < (1 - 2e)[(A - r)is - A)]i/2[1 + \{s + rf{M\Y^C^{J)\ (3.21) 

o 

for some positive constants Ci(J) and C2{J) which depend only on J. 

Note that the operator T has a spectral gap (r — A + 7^e(iV), s — A — 7^e(A^)) =: (— d+). 
Denote w{X) = (A — r)(s — A). Evoking the above formulas we obtain 

ci+ > s-A-Ci(J)w;(A)(inf A„)-^ (3.22) 

p>N 

and 

d- > X - r - Ci{J)w{X){mf Xp)-\ (3.23) 

p>N 

Let r{N) = (infp>7v Ap)~^. Inequalities (13.211) . (13.221) and (13.231) allow to compare and d+-d_. 
We find 

d+-d_-(3^ > {s - X - Ci{J)w{X)r{N)){X - r - Ci{J)w{X)r{N)) 
- (1 - 2e)Mm + lis + rfC,{J)r{Ny]' 

o 

= w{X) [(1 - Ci(J)(A - r)r(iV))(l - Ci( J)(s - A)r(iV)) 

-(1- 26)2(1 + i(. + r)2C,(J)r(iV)^)] 
> w{X)[Ae-Ae^-C,iJ)r{N)l 

for some positive constant C^lJ). 

Choosing so large that r{N) < 4e2(C3(J))~^ the above inequalities show that 

d+-d_-p^> 4M;(A)e(l - 2e) > ew(A), 

for e e (0,1/2). 

Therefore Lemma [3.31 implies that J — A + ^4(7) is invertible and 

||r^(J-A)-V|| = ||(J-A + A(7))-^|| 

< [A+-(A2^ + /32-rf^.rf_)V2]-i 

< 2A+[(rf+-rf_-/52)]-i 

< 2A+{ew{X))-^ = {s - r){w{X)e)-\ 

Using the last inequalities and repeating the reasoning given at the end of the proof of Theorem 
12.11 we get the desired estimate, thus completing the proof of Theorem] 



Remark 3.5 The constant rj in Theorem 12.11 can be made arbitrary small. In turn the choice 
ofr7 = l/2 — ein Theorem 12.21 is optimal as will be shown below in Example 14. 1[ 
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3.3 Proof of Theorem [2731 

Fix A such that 3? A < d and e G (0, 1). Choose such that 

infA„>l, 2(inf A„)-^exp[(l-e)2(rf-3?A)] < e. (3.24) 

n>N n>N 

Fix an arbitrary A^^i > A^. Now define the sequence p{n) by 

0, n< N, 

Pin)={ E::^(v^)-^ N<n<N„ 



Put 7 = (1 — e)y/d — 9fJA and consider the multiphcation operator (pNi = in £^(N). Note 
that (pN^ is bounded and invertible (for any Ni > N). By repeating the computation given in 
the proof of Theorem 12. II we have 



- J =: A(7) (3.25) 



is a finite rank Jacobi operator with entries and bk (see the proof of Theorem 12. ip . where 
ai = • ■ ■ = ttN-i = bi = ■ ■ ■ = &Ar-i = and = 6^ = for A; > A^i. By definition of ^4(7) one 
can easily check that 

||ReA(7)|| <72(l + 2e^'(inf Ap)-^) (3.26) 

p>N 

Define the operator A in on the domain D{J) {D{A) = D{J)) by 

Af := J/ + ReA(7)/ + ^ImA(7)/. 
Hence for complex such that 3fJ/i < d — ||Rey4(7)||, A — /i is invertible and 

II (A - /i)-l < [{d - ||ReA(7)||) - ^^^r' (3.27) 

Since 

A = J + A{j) = 0-;j0jVi 



for A fixed as above we have 



(A-/i)-i=0^i(J-/i)-V^,. (3.28) 



In turn (13.281) implies that 

K(A-A)-V,e„)| = |((J-A)-V7v,ei,0^^e„)| = 0^;(n)|((J-A)-V,e„)|, 

for any N < n < Ni. 

Using inequality (13.271) we obtain 
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provided that < n < iVi and 3?A < c/ - ||Re A(7)||. 

Finally, using the definition of 7 as well as fl3.24p and (13.261) one can readily check that 

||ReA(7)|| < (1 -e)(rf-3?A). (3.30) 

Combining the above relations (13.291) and (13.301) and the fact that A^i was arbitrary com- 
pletes the proof of Theorem 12. 3[ 

Remark 3.6 Note that the above = N{^X,e) tends to 00 as either 3?A goes to —00 or 
e — > 0. In the interesting region {A : 3fJA ~ d} the dependence of on 3fJA disappears. 



4 Optimality of the decay bounds 

In this section we discuss two specific models illustrating that the bounds obtained above are 
optimal. 

Example 4.1 The first model of unbounded Jacobi matrices J concerns Theorem 12.21 and is 
given by 

qn = 0, Xn = n + Cn, n G N, (4.1) 

where (c^) is a two-periodic sequence (ci, C2, ci, C2, . . .) such that ci 7^ C2 and A„ > for all n. 
Unbounded Jacobi matrices quite frequently fall into one of two extreme cases, namely that 
either aess{J) = M or cTessiJ) = 0- What makes the class (14. ip interesting is that it is non-trivial 
in this respect. It can be shown that (TessiJ) = {~\p\Ap\)^ where p = ci — C2. In fact, the 
spectrum of J is purely absolutely continuous in ]R\ [— This is proven in [8jand [T2] by 
finding asymptotics of solutions of the equation Ju = \u and using the method of subordinacy 
(for the a.c. spectrum) for this and related classes of entries. In the same papers it is shown 
that o"(J) in (— |p|, |p|) is empty. In other words we have exactly the situation considered in 
Theorem 12. 2[ with r = —\p\ and s = |p|. 

Moreover, in [8j and |T5] it was shown that for A G (r, s) there exists a solution Ui{n) of the 
equation 

Xn-iu{n — 1) + qnu{n) + Xnu{n + 1) = Au(n), n > 1 (4.2) 



such that 



ui{2n - 1) 
Ui{2n) 



[(J - A)"^ei](n) = aMi(n), n > 1, a G M (4.3) 

(4.4) 



j =d„exp[-v/;^^^^(2A;)-i]^(e_ + o(l)) 
^ fc=i 

where dn = (-1)" nLi(l - V2A;), e_ = ° ^ , and 



S :-- 



1 1 

w+r — ui+r 
A ' A 



Comparing (14.40 with the estimate of Theorem 12.21 completes the proof of the above men- 
tioned sharpness of the result. 
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Remark 4.2 The above solution ui is the unique £^-generahzed eigenfunction to which the 
bound found in Theorem 12.21 apphes. As dn ~ n~^/^ is subexponential, the decay of ui at +oo 
is governed by exp ^— a/p^ — X]fc=i(2^)^^ j • With r] = 2~^, the n-dependence of the exponent 
through the factor Ylk=i^^^ ~ J2k=i^/'^k and (for A close to the boundary of the essential 
spectrum) the dependence on the distance of A to the spectrum given through — = 
— A)(A + p) are exactly as found in Theorem 12.21 

The next example we present below illustrates that the estimates of Theorem 12.31 are also 
sharp. 

Example 4.3 Consider the sequences given by 

Qn = -2n, Xn = n, n e N. (4.5) 

Note that J with entries defined by (14. 5 p is bounded from above by — /. 

Fix A > — 1 and e G (0, 1). By the Birkhoff-Adams theorem [6] there exists a basis u±{n) of 
solutions of (14. ip with asymptotics given by 

u^{n) = n-^/Sxp(±2v/(A + l)n)(l + o(l)). (4.6) 

Applying Theorem 12.31 (in this case for J bounded from above by (— /)) there exists N such 
that for n > N we have 

n 

\u.{n)\ < (e(l + A))-^xp[-(l-e)v/rTA^(v^)-^]. 

k=N 

Since X]fc=Ar(V^)^^ — 2v^) by comparing (14.60 and the last estimate we conclude that the 
value 1 — e in the formula for 7 (see the the proof of Theorem 12. 3p is optimal. 

Problem. In case that J is a bounded Jacobi matrix it is well known (see Theorem 2.3 in [18]) 
that the spectrum of J as an operator in P does not depend on p. However, for unbounded 
J this result does not apply. Nevertheless, in the case yZ'^-i t- = 00 the estimates given in 
Theorem 12.21 imply that cr^(t/) ^ <^^{J) for any p G [l,C)o], where o"P(J) is the spectrum of J 
considered on the maximal domain in This can be easily seen by applying the Schur test to 
((J — A)~^ej, e„). Does the opposite inclusion also hold true? 



5 A criterion for the existence of square-integrable so- 
lutions 

While the results of this section could be stated and proven for general Jacobi matrices, we will 
for simplicity assume that J is given by (12. ip with zero- diagonal, = for all n G N. 

We will compare J with a second Jacobi matrix Jq, also with zero-diagonal and weights A°, 
n G N. Both weight sequences satisfy (12.30 to guarantee self-adjointness of J and Jq. 
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We assume that the weights of J and Jo coincide on a sequence of increasing intervals: Let 
Xfc G N, xi < 0:2 < < . . ., £fc G N with hnifc^oo^fc = 00 and Xk + ik < Xk+i — ik+i for all k, 
and assume that 

A„ = A° forallnGUfc{a;fc-4-2,Xfc + 4 + l}. (5.1) 
Denote 4 := {xk - ik, Xk + ik} and 

Afc := max{A„ : - 4 - 2 < n < + 4 + !}• 

Theorem 5.1 Suppose that for all 7 > 0, 

00 

J2 Afc(Afc-i + Afc + Afc+i)e-^^'=/^'= (xfc+i - Xfc_i) < 00, (5.2) 
fe=i 

i/ien /or almost every E E M \ cr(Jo) ^/iere exists a non-trivial square- summahle generalized 
eigenf unction of J to E. 

Note that (15. 2p allows for situations where the distances between centers of barriers Xk+i—Xk 
can grow significantly faster than their size 24- This leads to applications where J may contain 
a lot of spectrum outside the spectrum of Jo, including entire intervals. However, this spectrum 
can not be absolutely continuous and will typically be pure point, as follows from the following 
well-known general result. 

Proposition 5.2 Let J he an unbounded self-adjoint Jacobi matrix in and I an open 

subset of M. Assume that for almost every E & I there exists a non-trivial square- summahle 
generalized eigenfunction of J to E. Then 

(a) aac{J) n / = 0, 

(b) (Jci-Jx) n / = for almost every A G M, where Jx = J + A(-, ei)ei. 

A proof of (b) for the case of discrete Schrodinger operators (i.e. all An = 1), using spectral 
averaging over the coupling parameter A, is given in [TSj. Due to the limit-point condition 
(12. Sp this proof extends to our setting. Part (a) follows from (b) as the absolutely continuous 
spectrum is invariant under rank one perturbations. 

The remainder of this section will be devoted to the proof of Theorem 15.11 Note first that 
it suffices to show the existence of a non-trivial square-integrable generalized eigenfunction for 
almost every E G where the compact interval is disjoint from cr(Jo), as M \ o"(Jo) 

can be exhausted by countably many such intervals. 

For E G M, > and z := E + ir] let 

«^(n,E) := ((J-^)-^ei,e„), riGN (5.3) 
be the Weyl-solution, i.e. the unique 4-solution of 

G N (5.4) 
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satisfying the boundary condition M^(0,i?) = —1. The m-function is given by m{E + irj) : = 
M^(l, E) and it is known that m{E+iri) — > m{E) as ?7 — exists and is non-zero for almost every 
G M. For arbitrary 2; G C let ^J^-^z) and ^{-^z) be the solutions of (15.41) with %Ij{Q,z) = — 1, 
ipil, z) = 0, ip{0, z) = and ip{l, z) = 1. Then 

Ur,{n, E) = ip^ri, E + irj) + m{E + iri)ip{n, E + irj), n >0. 

We have pointwise in n and E that ijj{n,E + irj) 'ip{n, E) and if{n,E + it]) <f{n, E) as 
1] ^ 0. Thus, for almost every G M, 

u{n, E) := -^{n, E) + m{E)if{n, E) (5.5) 

exists and is a solution of (15.41) with z replaced by E. Theorem 15.11 is proven if we can show 
that, for a suitable tjq > 0, 

00 

sup |M^(n, i?)|^ < 00 for a.e. G [a, (5.6) 

o<v<vo 

as this implies that u{-,E) is square- summable. One has u{l,E) = m{E), thus u{-,E) is also 
non-trivial for almost every E. 

In the following we will find bounds for the £^-norm of separately within the barriers Ik 
and for the intervals between any two given barriers. 

To handle the barriers, define ik by rounding down ik/^ to the closest integer and let 
h '■= [xk — ^k,Xk + f-k]- Also, write Xk and Xk for the characteristic functions of Ik and Ik- A 
calculation using (15. ip shows that 

{-\n~iUn{n-l), n e {xk - ik,Xk + ik - I}, 
-XnUrf{n + l), ne {xk + ikjXk + ik + 'i-}, (5.7) 
0, else. 

We write 

XkUri = XkXkUrj = Xk{Jo -{E + iv)y^XuAJo - iE + iri))xkUn, (5.8) 

where 

Uk ■■= {xk - 4 - 2, Xfc - 4 - 1, Xk - 4, Xk + 4, Xk + ik + 1, Xk + ik + 2}. 
Plugging (15.71) into (15. 8p yields 

< 4A^||xfc(Jo-(i5 + ^r7))''xi/JI' 5^ |wr,Hr 

< 4A^a^ 5^ K(n)p, (5.9) 

where the elementary bound (|a| + |6| + |c| + \d\y < 4(|ap + |6p + |cp + |c/p) was used and we 
have set 

flfc := sup WxkiJo- {E + ir])y\uJ■ 
o<v<■no,Ee[a,b] 
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We now consider the intervals between barriers. For this let Jk be the finite Jacobi matrix 
found by restricting J to the interval [xk,Xk+i]- We have, also considered after restricting to 

[Xk, Xk+l], 

(1 - Xfc - Xk+i)ur, = (Jk - iE + ir])y^{Jk - {E + ir])){l - Xk - Xk+i)ur,. (5.10) 
Similar to (15.71) we find 

{-Xn-lUr,{n - 1), ne {Xk + ik + l,Xk + ik}, 
-XnUr^in +1), n G {Xk+l - 4+1 - 1, Xk+l - 4+l}, 
0, elsewhere in [xk,Xk+i]- 

(5.11) 

With 

Vk '■= {xk + 4 — 1, a^fc + 4, Xk+l — ^k+i,Xk+i — 4+1 + 1} 
we get from fl5.10l) that 

+ 1 — "^fe+l — 1 

E l^'^l' ^ 4(A^ + AUMJk -{E + ^v)r'r E l"^''!'- (5.12) 

n=Zfc+4+l "^^fc 

Given that C U Jfe+i we can bound the term ^^eVk l^*?!^ right hand side of ( 15. 12^ 

by two terms of the form (15.91) to arrive at 



El |2 ^ 16(A| + A^_^;^) / . 2 2 I |2 I A 2 2 | | 

'''''' AkiEV [^kakl^Kl +Vi«fc+i 2^ K\ 

n=Xfc+4+l \ '^eC/fc neUk+1 

Here we have also used that 

1 



(5.13) 



\\(J,-{E + z7^))-'\\ < 



ME) 



abbreviating Ak{E) := dist(£', cr( J^)). 

Ultimately, we can bound the £^-norm of over the entire interval [xk,Xk+i] by the sums 
in (15. 9p (for Ik and h+i) and in (I5.13P to get 



YlKinf < 4A|ai(l+il^^l||ili)5:i„,(„) 

(^l + ^^^^J^E K(n)l • (5.14) 
Summing over k > K for any given i^T G N yields 



I A2 

+ 4AUial^, (1 + 



ne(7fc+i 



n=a;if k=K n=X]^—if^—2 
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where 



,22/1 I ^^k-l^-^^k I 



WB):=Ai<.Ul + ^^ + ^^l. (5.16) 



We will now use 

Lemma 5.3 Under the conditions of Theorem \5.1\ it holds that 

lim hk{E) = for a.e. E e [a,p]. (5.17) 

Before proving Lemma 15. 3[ we show how it is used to complete the proof of Theorem 15. 1[ Let 
E G [a, j3] be such that bk{E) — > and also m{E + irf) m{E) G C as 77 — > 0. By Lemma [5^3] 
and our remarks at the beginning of the proof, this holds for almost every E G [q;,/3]. 

There exists ko{E) such that bk{E) < 1/32 for all k > koi^E). By our general assumptions, 
any given n G N is not covered by more than two of the intervals [x^ — £fc — 2, + £fc + 2]. 
Thus, if we pick K = kQ{E) in 05.151) . it follows that 

tiGN n=l n=XK 

XK-l ^ 

n=l neZ 

or 

XK-l 



$^K(n)r<2 5^ \u,{n)\'. (5.18) 



neN n=l 

The crucial fact is that this holds uniformly in rj. By (15.51) . limM^(n) = u{n) exists for each of 
the finitely many n G [l,x;<: — 1]. Thus the right hand side of (15.181) is uniformly bounded in 
Tj G (0, 1]. This concludes the proof of (15. 6p and thus of Theorem 15.11 

We finally prove Lemma \5.3[ where Theorem 12.41 is the crucial ingredient. The interval 
[a,P] has positive distance do := dist([a, /?], (t( Jo)) from the spectrum of Jq. At this point 
we choose rjo := do/ 8, meaning that tjq < a/(A — r)(s — A)/8 for all A G [a,/?], where r : = 
max(cr(Jo) fl (—00, a)) and s := min(cr(Jo) fl (/9, 00)). Thus, by Theorem 12.41 there exist 
constants C < 00 and 70 > such that 

a;fc-4-l \ / ^k+^k ^ 

ak<C exp I -70 Y Y] ^'^^^^\ ^ A~ 

n=a::fc-4-l J ^ n=Xk+ik " 

Here we have split the four points 'va.Uk into the two pairs to the right and left of and applied 
(12. 8p separately. Using the definition of and that Ik < we arrive at 



(5.19) 
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Let Ak := {E G [a,P] : Ak{E) < ak} with ak to be chosen later. has Xk+i — Xk + 1 
eigenvalues. Thus \Ak\ < 2ak{xk+i — Xk + I) < 4:ak{xk+i — Xk). Suppose that 

^ak{xk+i - Xk) < oo. (5.20) 

k 

Then \ Ak\ < oo and it follows from the Borel-Cantelli lemma that for almost every E G [a, (3] 
there exists kQ{E) G N such that Ak{E) > ak for k > ko{E). For such E and k we conclude 
from fl5:T6D and flCTD that 

/ A2 I /\2 A2 I A2 \ 

bk{E) < CA^e-^°^^/2^^ ( 1 + ' + ' ) . (5.21) 

Let 7i G (0,7o/2) and 

«^ := A^(A^. + ALi)e-^^''=/^'= + Ati(A^. + Ati)e-^^^-^/^'=-. 
With this choice of ak, (15. 2p shows the summability of \Ak\. Moreover, from (I5.2ip we get 

bkiE) < C (^A^e-^o^*/'^^ + e-(37o-7i)4/A, j ^ 

which tends to as A; — > oo by (15.21) . Lemma [5.31 is proven. 

6 A class of Jacobi matrices with a mobility edge 

As an application of Theorem 15.11 we now provide a class of unbounded Jacobi matrices, which 
exhibit a transition from spectral regions with purely absolutely continuous spectrum to a region 
with dense pure point spectrum. Due to the connection of spectral and transport properties 
this is called a mobility edge in the physics literature. Our example is motivated by classes of 
bounded Jacobi matrices considered in [21] which exhibit the same behavior. Here, as well as 
in [21], a purely absolutely continuous Jacobi matrix is subjected to a slowly oscillating pertur- 
bation, which generates regions of dense pure point spectrum. As opposed to the examples in 
|2T] , in our example an additional periodic modulation sequence is needed to open up a gap in 
the purely absolutely continuous spectrum of the unperturbed unbounded Jacobi matrix with 
weights n", < a < 1. 
Define the weights by 

A„ = n" + c„(/^(n^). (6.1) 

Here : M ^ M is twice continuously differentiable and periodic, i.e. (p{x + T) = ip{x) for some 
T > and all a; G M. Furthermore < ip{x) < 1 for all x, inf = and sup(/? = 1. The 
sequence c„, n G N, is 2-periodic, c„ = Ci for all odd n and c„ = C2 for all even n with ci > 0, 
C2 > and Ci 7^ C2. Write c := |ci — C2I. 

Theorem 6.1 Let J be the Jacobi matrix on £^(M) with zero- diagonal and weights given by 
h6. Suppose that < a < 1 and < 7 < (1 — a)/2. Then 
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(a) a{J) = R, 

(b) J is purely absolutely continuous m R \ [— c, c], 

(c) J has no absolutely continuous spectrum in (— c, c), 

(d) J\ := J + A(-, ei)ei has pure point spectrum in (— c, c) for almost every A G M. 

We start with the proof of part (a). Denote by Jq the Jacobi operator with zero diagonal 
and the weights An^"* = n", a G (0, 1]. Let {Jou)n = Xun, n > 1, A G M. Using Theorem 3.2 of 
[TT] we know that 

\un\ ^ n~°'^'^ as n oo. (6.2) 



Let Ei := i ~' , where < 6 < Assume that (p{xo) = and consider the sequence 

{xi} C M such that x] = Xq + iT. It follows that Xi ~ T^hj}h ^ 

For rij := [xj] choose a sequence {Aj} C N of even numbers such that 



A, > e,n]-^ (6.3) 

and 

Ai < Mn]'^~' (6.4) 

for some M > and < £ < 1 — a — 7 — 7(5. Consider the sequence of intervals Jj := [nj, nj + Aj]. 
Denote fij := rii + Aj/2. We let /?j := 2/Aj and claim that the sequence {ti^*^} G £^ given by 

[0, n^h 

[ Un[l - (3i{n - hi)], hi <n <ni + Aj, 

is a Weyl sequence for J at the point A: 
By definition of f one can verify that 

||i;«||2~(A,nr")V2. (6.5) 
Moreover, observe that the above choice of Aj implies 

sup |v?('^''')| —^0 as 2 — > 00. 

In fact 

\^{n^)\ = \v{n^)-v{x])\<y\U{n, + A,y-n]\ 
< Cn]—< CMn-' 

m 

for some C > and i sufficiently large. Here we have used (16. 4p . Hence 
for some Mi > and large i. 
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We find for Ui < n < Ui (tlie "+" -signs in tlie following calculation) and for rii < n < n-j + Aj 
(the "—"-signs), respectively, 

[(Jo - X)v^\n) = [±{n - l)°M„_i(ra - 1 - Ui) ± n"M„+i(ra + 1 - Ui) - X{±u„){n - hi)]l3i 
= ±Pi{n - hi)[{n - l)°u„_i + n°'Un+i - Am„] 

±[{n - l)°'/3i{-Un-l) + n°'Un+l/3i - X/3iUn]. 

Note that the first term in the last expression vanishes as JqU — Xu = 0. The only other values 
of n such that [( Jq — A)f '^*^] (n) ^ 0, n = rii, hi, ni + Ai, give slightly different expressions which 
do not contribute significantly to ||(Jo — A)u^*^||. All this and the bound fl6.2p lead to 

II (Jo - A)t;«f < M^irii + A,)2-/52,,--(A, + 3). (6.7) 

Combining ( 16. Sp . ( 16. 6p and ( 16. 7p we find 

||(J- A)t;«ll 



for some M3 > and large i. From ([631) it follows that ^nf < M4~^ as i ^ 00. This 
completes the proof of part (a). 

(b) Fix A G {x G M : |x| > c}. We want to study the behavior of generalized eigenfunctions 
of Ju = Xu. Denote by Bn the transfer matrix of J given by 



1 

An-l A 



An Xn 



As in our previous works [ini[ITl[l2] it is of advantage to compute the products i?2n-B2n-i- We 
have 

A2n ^ ^ l + C2V5((2n)7)(2n)- ^ ' 

^ . 1 + (2.)---^((^- - f - - ^);) - --(^-)°" + (6.9) 

A2n-i ^ ^ (l-l/2n)" + ci</.((2n-l)^)(2n)-" ^ ' 

where p„, 2;„ G l^. Since < 7 < (1 — a)/2 and ip' and ip" are bounded, the sequences 
{¥?((2n-2)T)}, {v?((2n- 1)^)} and {(/?((2n)T)} belong to the class D^{{2n)"), see Example 2.1 
in [H]. Recall here that for a non-negative weight sequence /i = the class -D^(/i) is defined 
in [9] as {x G : Ax G Z^(/i),A^x G /""^(yu)}, where a sequence x = {x„} is in if 
|a;(n)|^/i(n) < 00 and A is the forward difference operator, i.e. (Ax)„ = x„+i — x„. It is 
clear that {n~"} is also in D'^{{2n)°'). 

The class D^((2n)" is closed under multiplication and division (by a sequence separated 
from zero), see Lemma 2.2 in [H]. Therefore we can write (16.81) and (16. 9p as 

^ = 1 + (2n)-"[(ci - C2M{2n - ly) + x„] + (6.10) 
^ = 1 + (2n)""[(c2 - cM{2n - 2)^) + + s„, (6.11) 
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where {xn}, {Vn} belong to D^((2n)°), \imxn = lim?/, 
By the same arguments 

= (2n)-X, 

where {t„}, {wn} G D^((2n)"), hmt„ = hmw„ = 1. 
Using (leTTUD . (16TT|) and (1^321) we obtain 

B2nB2n-i = "/ + (2n)-"r(n) + 

with G The matrix V{n) has entries 

V^iiH = -{C2 - CiM{2n - 2y) - yn, 

V2i{n) = -Xwn{l + i2n)-"[{c2~CiM2n-2y + yn]}, 
V22in) = -[(ci-C2)v5((2n-l)^) + x„] + A2(2n)~Xti^„. 

Thus V{n) is in D^{{2n)'^). 
Note that 

discrV"(n) := {tTV{n)f - AdetV{n) 

= -A[X' + {C2 - ci) V((2n - iyM{2ny)] + 0(1), 

as n — s> oo. It follows that limsup discry(n) < 0. Thus we have verified all the assumptions 
to apply Theorem 5.1 of [H], and in particular the asymptotic formula (5.17) there. This shows 
that Ju = Xu has no subordinated solutions in the sense of [H]. It follows from the results of 
[H] that J is purely absolutely continuous in M \ [— c, c]. 

To prove parts (c) and (d), let Jo be the Jacobi matrix on £^(N) with weights = n"+c„. It 
was shown in [12] that Jo has at most finitely many spectral points in (— c, c). Thus (— c, c)\a{Jo) 
is the union of finitely many open intervals Uj. 

Fix a Uj and a compact interval [a,l3] C Uj. We will show that J has a non-trivial square- 
summable generalized eigenfunction for almost every A e [a,/?]. This implies that the assump- 
tions of Proposition 15.21 hold with / = (— c, c) as (— c, c) \ cr(Jo) can be exhausted by countably 
many such intervals. Thus (c) and (d) follow. 

Choose Xq G [0,T) such that ipi^Xo) = 1 and let 6 := dist([Q;, /?], ±c). There exists e > 
such that ip{x) > 1 — 5/(2 max{ci, C2}) for all x G [xq — e,xo + e]. For every G N, let Xk be 
the integer closest to the center of the interval [{xq + kT — eY^"', {xq + kT + eY^'^] and ik an 
integer approximately equal to a quarter of the length of this interval. Let J^ be the Jacobi 
matrix with weights 

_ f A„, n G Ufcfxfc - 4 - 2,Xfc + 4 + 1], 

' ri" + c„, else. 

Then |A^ — (n° + c„)| < 6/2 for all n and thus [a,P] is contained in a spectral gap of J^. We 
will complete our proof by showing that the assumptions of Theorem 15.11 are fulfilled with J^ 
in place of Jo. 



0, and {r„}, G 



(2 



n) 



Wr. 



(6.12) 
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For this, note that Xk ~ {ko + kTy/^ and that Xk+i — ~ "^{xq + kTy^"' ^ and ~ 
{xq + kT)°'^"' are polynomially bounded in k. On the other hand, ik ~ ^{xo + kTY^"'^^ and 
thus 

4/Afc (xo + A;T)'t'-\ 

7 

The latter is polynomially growing in k since 1 — a > (1 — a)/2>7. This implies that (15.21] 
holds for all 7] > 0. 
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